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Loop quantum cosmology is an application of recent developments for a non-perturbative and background 
independent quantization of gravity to a cosmological setting. Characteristic properties of the quantization 
such as discreteness of spatial geometry entail physical consequences for the structure of classical singular- 
ities as well as the evolution of the very early universe. While the singularity issue in general requires one 
to use difference equations for a wave function of the universe, phenomenological scenarios for the evolu- 
tion are based on effective equations implementing the main quantum modifications. These equations show 
generic bounces as well as inflation in diverse models, which have been combined to more complicated 
scenarios. 



1 Introduction 

The universe is, on large scales, well described by general relativity which provides the basis for mathemat- 
ical models of the possible behavior of a universe. From solutions of Einstein's field equations one obtains 
the geometry of space-time once the matter content has been specified. On this classical level, however, 
the description will always remain incomplete as a consequence of singularity theorems: any space-time 
evolved backward in time from the conditions we perceive now will reach a boundary in a finite amount 
of proper time. At this boundary, the equations of the theory, and thus classical physics, break down, often 
accompanied by curvature divergence. For the development of complete universe scenarios the classical 
theory of general relativity thus has to be extended. 

Such an extension is often expected to come from quantum gravity, where not only general relativistic 
but also quantum effects are taken into account. One approach, which is background independent and 
non-perturbative and can deal with those extreme conditions realized at classical singularities, is loop 
quantum gravity (Q. It is a canonical quantization and turns the classical metric and extrinsic curvature 
into operators on a Hilbert space. The classical geometrical structures are thus replaced by properties of 
operators which by itself leads to a different formulation. While classical geometry has to re-emerge as an 
approximation on large scales, on small scales quantum behavior has to be taken into account fully. (For 
general aspects of quantum theory in the context of cosmology see |2|.) Then also the singularity problem 
appears in a different light as the basic object is not a space-time metric which cannot be extended beyond 
singularities in the classical evolution but a wave function. The wave function, also, is subject to equations 
which need to be analyzed in order to see if it always gives a complete solution telling us what happens at 
and beyond classical singularities. 
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Classical space-time is thus replaced by quantum space-time whose effects are most important on small 
scales such as those of a small universe close to a classical singularity. A fully quantized system is indeed 
necessary to describe states right at a classical singularity, but this is in general complicated at technical 
as well as conceptual levels. Close to classical singularities it is thus helpful to have effective systems 
which are of classical type, removing interpretational issues of quantum theories, but take into account 
some quantum effects. With those ingredients it is then possible to develop several complete scenarios for 
universes, and at the same time obtain potentially observable phenomenological effects. 

2 Loop quantum cosmology 

Just as Wheeler-DeWitt models, loop quantum gravity is a canonical quantization, i.e. it starts with a fo- 
liation of space-time into spatial slices S t . Canonical variables are the spatial metric q a f, and its momenta 
related to extrinsic curvature K a b of the slice Q. The canonical Poisson algebra then is to be turned into 
a suitable operator algebra and represented on a Hilbert space. As always in field theories, however, oper- 
ators for the field values in single points do not exist and the fields have to be smeared first by integrating 
them over extended regions. For field theories other than gravity, this is usually done on 3-dimensional 
regions using the background metric to define an integration measure. But for gravity, the metric itself is 
dynamical and thus to be smeared, and using an additional metric for smearing would introduce a back- 
ground. A different procedure is required for gravity, i.e. we need to smear fields but do so in a background 
independent manner. 

2. 1 Holonomy-flux algebra 

Indeed, the success of a quantization procedure often depends on the choice of basic variables. For gravity, 
it is most helpful to transform to new basic fields, given by Ashtekar variables |4j[3]. These variables are 
also canonical, but instead of the spatial metric one uses the densitized triad Ef, related to the spatial metric 
by EfE° = dct qq ab , and the Ashtekar connection A l a — T l a + jK^ defined with the spin connection F l a 
and extrinsic curvature components K l a . In addition, there is the Barbero-Immirzi parameter 7 > Pl|6). 
which for simplicity will be set equal to one in what follows. (Its value can be computed from black hole 
entropy and is smaller than but of the order of one 1 7 1 .) In those variables, the extrinsic curvature term in A l a 
makes it canonically conjugate to the triad, and T l a provides the transformation properties of a connection. 
For Ef, the important properties compared to the metric are that it also knows about the orientation of space 
since it can be left- or right-handed, and that it is dual to a 2-form eabcEf. The orientation will be essential 
later on in the discussion of singularities, while the transformation properties of a connection and dual 
2-form, respectively, are important right now because they allow a natural smearing without introducing a 
background metric. 

We can then integrate a connection along curves, where for good gauge properties we also take the path 
ordered exponential, i.e. use holonomies 

h e {A) = Vexs> J n A? a e a dt (1) 

for curves e C X and fluxes 

F S (E) = f T*E?e abc dy a dy b (2) 
Js 

for surfaces S C S. (SU(2)-generators Tj = —\i<Jj with Pauli matrices <jj appear because this is the 
gauge group of triad-rotations not changing the metric.) There are then no 3-dimensional smearings, but a 
1- and a 2-dimensional one, which, as it turns out, adds up to the right overall smearing for a well-defined 
quantization. 
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This is the basis of the background independent quantization provided by loop quantum gravity |8|, 
and it has many crucial properties as direct consequences. First, the holonomy-flux algebra, defined by the 
Poisson relations, is, under weak mathematical conditions, represented uniquely on a Hilbert space together 
with a unitary action of diffeomorphisms of E 0. The latter property is required for the independence of 
gravity under the choice of spatial coordinates. This representation is cyclic, i.e. there is a basic state from 
which all others can be obtained by repeatedly acting with operators of the algebra. Other characteristic 
properties of the basic holonomy and loop operators are then very different from those in a Wheleer- 
DeWitt quantization: There is no operator for connection components or even their integrals, but only 
for holonomies unlike in a Wheeler-DeWitt quantization where extrinsic curvature components are basic 
operators. Moreover, flux operators have discrete spectra and so do geometrical operators such as area 
and volume 1101 . With these properties, one can then use the basic representation to construct classes of 
well-defined Hamiltonian constraint 1 1 1 |and matter Hamiltonian operators 1 12 1. 



2.2 Isotropic quantum cosmology 

These techniques also make a symmetry reduction possible which is much closer to the full theory than 
a Wheeler-DeWitt model would be 1 13]. When symmetries are imposed, natural sub-algebras of the full 
holonomy-flux algebra are defined which, using cyclicity of the representation, induce the basic represen- 
tation of a model. Here, one simply acts only with those operators in the distinguished sub-algebra for 
a given symmetry, and thus obtains less states than one would get from the full holonomy-flux algebra. 
In this sense, the basic representation of models, which is so important for other physical properties, is 
directly obtained from the full quantum theory: quantization is done before performing the symmetry re- 
duction, at least as far as the basic representation is concerned. More complicated operators such as the 
Hamiltonian constraint can then be constructed from the basic ones following the steps in the full theory 
by analogy. 

The simplest case, isotropy, serves as a good example to illustrate the basic properties 1 14 1. Classically, 
there is a single gravitational degree of freedom, the scale factor a with its momentum p a = — 3(47rG) _1 ad 
(G being the gravitational constant). These are components of the spatial metric and extrinsic curvature, 
which now are replaced by Ashtekar variables. In the isotropic case, there is again a single canonical pair 
(c,p) with 

\p\ = a 2 , sgn(p): orientation (3) 

and 



c=i(fc + d) , k = 0: flat, k = 1: closed. (4) 

These variables could be quantized directly in a Wheeler-DeWitt manner, i.e. a as a multiplication 
operator and p a — —ihd/da on square integrable functions (of only positive a, which means that p a in this 
manner is not self-adjpoint), but an analog of the latter operator does not exist in the full theory. Through 
the induced holonomy-flux representation it is rather the exponentials exp(i/ic/2) for any real fi (to be 
thought of as related to the parameter length of a curve) which are basic in addition to the densitized triad 
component p proportional to a flux. The induced representation is then defined on a Hilbert space with 
orthonormal basis of states 

( C | M ) = e^ c/2 , fx e K (5) 

and basic operators 

p\n) = ±4,i| M ) (6) 

e5V2|p) = Im + m')- (7) 
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Note that p is self-adjoint and e^' c / 2 is unitary since the full range of real numbers for /j, is allowed taking 
into account the orientation freedom in a triad. In this representation, basic operators of the model have the 
same properties as those in the full theory 1 15 1: p has normalizable eigenstates and thus a discrete spectrum 
(nonetheless, the set of eigenvalues is the full real line, which is not in conflict with the discreteness of the 
spectrum because the Hilbert space is non-separable) and one can see that there is no operator for c but 
only for exponentials e vc / 2 not being continuous in /x. 

Compared to a Wheeler-DeWitt quantization, however, the properties are very different. In that case, 
the scale factor, related to p, would have a continuous spectrum and extrinsic curvature p a , related to c, 
would be represented directly as an operator. The properties realized in loop quantum cosmology as in the 
full theory are a consequence of strong restrictions coming from background independence and its transfer 
to symmetric models. 

On this basic representation we can construct more complicated operators, most importantly the Hamil- 
tonian constraint. Since the properties of basic operators are as in the full theory, the construction can be 
done in an analogous manner, only adapting to the symmetric context where needed. Properties of com- 
posite operators are then also close to those in the full theory, although for them the relation, as of now, is 
not as tight and there is no derivation of symmetric composite operators from the full ones. 

For the Hamiltonian constraint, we start from the classical Friedmann equation 

H = -6 [2c(c - k) + k 2 ] y/\p\ + 8wGH mme r(j>, <f>,Pt) = 
with matter Hamiltonian -ffmattcr, e.g. 

= \a- 3 pl + a 3 V((f>) (8) 

for a scalar with momentum p^ and potential V(<p), written down in isotropic Ashtekar variables. This 
is, using the volume operator V = \p\ 3 ^ 2 with eigenvalues 

= (4ImI/6) 3/2 (9) 

and replacing factors of c by exponentials 0, quantized to an operator equation H\ip) = for states 
\tp) = Y^n ^vWIm) given by a difference equation H 1 41 1 1 61 IT31 such as 

(V^ +5 - ^ +3 )e 4fc W+4(0) - (2 + fc 2 )(V^ +1 - ViWcW (10) 

+ 0V3 - ^-5)e- lfc VV-4(^) - -f^G4ff ma ttcr(Ai)^W>) 

for VvW- Unlike Wheeler-DeWitt equations which are differential, there is thus a difference equation 
as a result of discrete quantum geometry. Nevertheless, for large /i ^> 1 the difference operators can be 
expanded in a Taylor series provided that the wave function is sufficiently differentiable 1171 1151 . This 
property can be included in semiclassicality conditions 1 18 19 1, but in more complicated models the exis- 
tence of sufficiently differentiable solutions may not be guaranteed 1 20 1 . 

2.3 Classical singularities 

With the propagation equation for the wave function on minisuperspace we can now address the issue of 
singularities in quantum cosmology. Compared to the classical situation, the setup has changed since we 
do not have evolution equations for the metric in coordinate time, but an equation for the wave function 
on minisuperspace. In semiclassical regimes, the metric can be reconstructed from the wave function, 
for instance making use of observables (as completed for a free, massless scalar in [21 1). But it is not 
guaranteed that a classical geometric picture is suitable everywhere for the behavior of a universe. From 
the point of view of quantum gravity, the propagation equation for the wave function is more fundamental, 
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and so also the singularity issue is to be addressed at this level. The question then arises whether or not 
initial values in one classical regime of large volume are sufficient to determine the quantum solution on 
all of (mini)superspace including regions which can be interpreted as being beyond a classical singularity. 
If this is the case, quantum space-time would not be incomplete and thus non-singular. 

For isotropic quantum cosmology, this question can be answered immediately. Minisuperspace is, first 
of all, enlarged compared to the usual metric space in that we have two regions differing by their orienta- 
tion sgn(p) and separated by degenerate geometries. Classically, they are thus separated by singularites. 
In quantum theory, the configuration space for wave functions also knowns about orientation via 

sgn(fi), and a classical singularity would occur at /i = 0. Unlike the classical evolution, the quantum 
equation for the wave function uniquely yields the wave function on one side if we pose initial values on 
the other side. In this way, aspects of quantum gravity give us, first, a new region of minisuperspace and, 
second, a unique extension between the two regions. In this manner, quantum gravitational models are 
singularity-free \22\. So far, explicit constructions include anisotropic (Bianchi class A) models 1 231 1241 
and as inhomogeneous cases spherical symmetry and polarized cylindrical gravitational waves B25II26I . 

The criterion of extendability used here is the most direct and most general one for the singularity issue, 
both at the classical as well as quantum level. In quantum cosmology, other criteria have been used such 
as the finiteness of the wave function at a classical singularity. Also this is realized automatically in loop 
quantum cosmology because, thanks to discreteness, there is always only a finite number of computational 
steps between initial values and a classical singularity. Infinities in the wave function could then only arise 
from diverging matter Hamiltonians, which does not occur in loop quantizations as we will explain more 
explicitly in what follows. In this context, we will also discuss curvature divergence which sometimes 
serves as a singularity criterion. 

At a more intuitive level, we obtain the picture of a universe which starts in a collapsing branch and 
evolves through a phase where continuous geometry fails but discrete quantum geometry remains mean- 
ingful. At the transition, the orientation of space changes implying that the universe turns its inside out. 
What exactly happens during the transition depends on the concrete matter model. If one has a parity vio- 
lating matter Hamiltonian such as that of the standard model, there are changes between the two branches 
even in the equations of motion. Otherwise, the two sides generically are still different from each other 
depending on the initial conditions for the wave function. Even if one starts close to a classical geometry on 
one side, it may then happen that after evolving through a violent quantum regime a new classical geometry 
will not be recovered. This does, however, not happen in matter models with a free, massless scalar which 
can be treated completely 1211 . In this case, solutions which are semiclassical for large volume bounce at 
small volume and become semiclassical afterwards. For those solutions, there is no difference in behavior 
when the orientation is changed. 

As already noted, the difference equation can be well approximated by the usual Wheeler-DeWitt equa- 
tion on large scales. On small scales, the equations differ considerably, but in some cases one can still 
compare the effects of initial conditions. For the Wheeler-DeWitt equation, initial conditions have origi- 
nally been imposed at a = with the intention of removing the classical singularity by requiring the wave 
function to vanish there 1271 1281 . However, this turns out to be ill-posed as an initial value problem, i.e. in 
most models only the trivial solution of a vanishing wave function exists. In loop quantum cosmology, on 
the other hand, dynamical initial conditions have been derived from the difference equation 1 18 29 1 which 
in many cases are comparable to DeWitt's condition ?/>(0) = but making it well-posed in the discrete 
setting |30). For a closed model, one can see 1311 that the dynamical initial conditions are closer to the 
no-boundary proposal |32| than to the tunneling proposal 1331 . But also in the discrete setting, the situa- 
tion is more complicated in less symmetric models such as anisotropic |20| or inhomogeneous ones |26| 
where analogous mechanisms are more difficult to realize. The possibility of conditions following from 
the constraint equation also depends on the ordering of the Hamiltonian constraint operator which may 
change coefficients in the difference equation. The ordering for d 1 01 as also used in 1 18 1 is not symmetric, 
and using a symmetric one removes additional conditions for wave functions in isotropic models. Sym- 
metric orderings are often required for semiclassical issues or for computing the physical inner product, 
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and moreover for a non-singular evolution in inhomogeneous models 1261 . In such a situation, isotropic 
models loose their dynamical initial conditions, but some conditions remain in anisotropic models [34| as 
well as the inhomogeneous case. Here, however, the analysis of implications for the solution space is still 
very incomplete. A general mechanism to provide initial conditions for the wave function of a universe is 
thus still outstanding. 



2.4 Matter Hamiltonian 

For a complete cosmological model, we also need to know its matter content and quantize its Hamiltonian 
for quantum cosmology. This, now, does not only include the matter fields but also geometrical factors in 
a matter Hamiltonian which need to be turned into operarors. For instance for a scalar, the Hamiltonian 
H t jj = \a~ 3 p 2 ^ + a 3 V((f)) has to become an operator in the field values but also in the scale factor a. We 
thus need a quantization of the inverse of a, or p, in order to quantize the kinetic term. 

At this point, properties of the basic representation become important: As we have seen, a loop quan- 
tization leads to an operator p which has a discrete spectrum containing zero, and such an operator does 
not have a densely defined inverse. This seems to be a severe obstacle, but it turns out that well-defined 
quantizations do exist [35 1, just as well-defined matter Hamiltonians exist in the full theory 1 12 1. In quan- 
tizations, the most obvious procedure is not always the successful one, and also here one has to start from 
alternative expressions for a~ 3 which are identical classically but lead to different quantizations, 
as, e.g., 



3 \ 3/(2-20 

— p^inhjihjWp] 1 })] (11) 



a' 3 



8irGlj(j + l)(2j 



with parameters < i < 1 and j G ^N, using only positive powers of p and "holonomies" hi = e cr ' of the 
connection component c. Inserting the basic operators and turning the Poisson bracket into a commutator, 
this can directly be quantized to a well-defined operator with eigenvalues I 36ll37l 



3/(2-21) 

^-( jW^Win t^l ■ (12) 



As one can see, the eigenvalues do depend on the parameters j and /, unlike the classical expression. 
Rewriting in the above manner thus introduces ambiguities as it is expected for the quantization of any non- 
basic operator. Important properties are, however, robust. For instance, for any choice of the parameters 
we obtain the classical behavior of a~ 3 at large values p 3> a peak around /i* = j£j> and decreasing 
behavior on small scales reaching exactly zero for p, = 0. 

For larger j, the sum in the eigenvalues contains many terms, and it can be approximated by viewing it 
as a Riemann sum of an integral. In this way, we obtain the effective density 

d(a)^ := %)^ } 2) = a-^(3a 2 /,4) 3/(2 - 20 (13) 

with p(p) = 6p/£p and 

«(9) = |9 1 - i (lT2(('7 + l) i+2 -|'Z-ir 2 ) 

- j^q ((q + 1Y+ 1 - Sgn(q - l)\q - l^ 1 )) . (14) 

The approximation becomes better for larger j as the number of terms in the Riemann sum increases, which 
can also be seen in Fig.^ 
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Fig. 1 Discrete sets of eigenvalues of d(a) for two values of j and 1 = 3/4 compared to the classical 
behavior a" 3 , the approximations by d(a)~L' and small-a power-law approximations. 

Analogous constructions of d(a) exist, e.g. d(ai, 0.2,0,3) in anisotropic models with diagonal metric 
components gjj = of. In contrast to an isotropic context, these functions are not required to be bounded 
for all configurations, i.e. all aj (see, e.g., |38|). The isotropic situation is very special in the way the 
approach to the classical singularity at a = happens for which there is only one possible trajectory in 
minisuperspace (which can only be followed with different rates in coordinate time t). In less symmet- 
ric models, there is much more freedom in the approach such as aj (t) in anisotropic models, and not all 
possible configurations are realized along the generic approach. Also in the full theory, degenerate config- 
urations exist on which inverse volume operators have unbounded expectation values |39|. At this point, 
dynamical information has to be used to find out if curvature remains bounded along effective trajectories 
of universe models. 




Fig. 2 Curvature potential on the minisuperspace of a diagonal Bianchi IX model with one metric com- 
ponent held fixed. Curvature is unbounded on minisuperspace, in particular at large anisotropies. Close to 
classical singularities (center and diagonal lines at x = or y — 0), however, it remains bounded. 

This can be done in homogeneous models such as the Bianchi IX model 1 24 1 where curvature is un- 
bounded on all of minisuperspace (Fig.|2ji but bounded along the effective and even quantum evolution as 
given by the difference equation. On smaller scales, also deviations from the classical approach happen 
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which imply non-chaotic behavior as the curvature walls start to breakdown from quantum effects I40l l41l . 
From this, one can draw conclusions for the generic inhomogeneous approach and structure formation us- 
ing the BKL picture [42 1, but for definitive conclusions inhomogeneous models have to be studied in more 
detail and in particular perturbative schemes for structure formation. 



2.5 Effective equations 

Difference equations, in particular partial ones in the presence of matter fields or ansiotropies and inho- 
mogeneities, are difficult to deal with, and also interpretational issues of the wave function arise at this 
level. While the full quantum setting is needed to discuss singularities and the discrete evolution in their 
vicinity, on larger scales it is more direct to use effective equations of classical type. Those equations 
are differential equations in coordinate time, but they are amended by modifications to capture quantum 
effects. This is analogous to effective action techniques, and indeed there is a general scheme for effective 
classical approximations of quantum systems 1431 M4 1 which can be applied to quantum cosmology in a 
canonical formulation and reproduces the usual effective action results 1 45 1 obtained by expanding around 
free field theories or the harmonic oscillator |44|. 

This scheme is based on a geometrical formulation of quantum mechanics [46 1 where the Hilbert space 
is interpreted as an infinite dimensional vector space with additional structure. The infinite dimensionality, 
even for a mechanical system, is the crucial difference between quantum and classical physics. There are 
thus additional quantum degrees of freedom which in some regimes have to be taken into account (often 
appearing as higher derivative terms in effective actions). From the inner product of the Hilbert space 
one obtains a symplectic structure as well as a metric on the vector space such that it becomes Kahler. In 
addition, the quantum Hamiltonian defines a flow on the Kahler manifold. With the symplectic structure 
and the flow one has a canonical system equivalent to the quantum system and in general involving all 
infinitely many degrees of freedom. In some cases, however, it is possible to approximate the flow by a 
dynamical system on a finite dimensional subspace of the full space, giving rise to an effective classical 
system. This may involve only the classical degrees of freedom, but with correction terms in the evolution, 
or also additional ones related to, e.g., the spread of wave functions. The metric on the Kahler space is 
needed only for the measurement process in quantum mechanics and issues such as the collapse or overlap 
of wave functions. Since there is no external observer in quantum cosmology and only one wave function, 
the metric may be dropped which means that one could, as mentioned in 0, weaken the Hilbert space 
structure required usually. Indeed, the geometrical formulation of quantum mechanics provides a unified 
scheme for generalizations of quantum mechanics [46 1. 

In loop quantum cosmology, several related methods have been applied in order to derive effective 
terms for equations of motion, although a complete derivation is still unfinished. The geometrical scheme 
is developed for quantum cosmology in 1431 with an asymptotic series of correction terms for isotropic 
models. Leading orders of such terms have also been derived in 1471 and with WKB techniques in |48|. 
Some of these terms have been used in applications 1491 1501 15 II 1471 . but many effects already show up 
at the level where additional degrees of freedom are considered only to the lowest order. The order of 
differential equations is then the same as classically, but terms in the Hamiltonian do change. In particular, 
the quantum Hamiltonian requires, in the presence of matter or other curvature terms, a quantization of 
a -3 which is modified at the quantum level. Effectively, a~ 3 in the matter Hamiltonian is replaced by d(a) 



from O 152 1 : 



H <f ,(a) = -d(a)pl + a a V^) 



(15) 
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From the matter Hamiltonian we obtain Hamiltonian equations of motion for the matter field, which 
now change on small scales [ 16 1. For a scalar, we have the effective Klein-Gordon equation 

l = ta^^-a 3 d(a)V'(® (16) 

where the usual friction term —3cj)d/a changes its form. Back-reaction from matter on geometry is encoded 
in the Friedmann equation which also changes. Substituting the effective matter Hamiltonian, we obtain 
the effective Friedmann equation 

a(d 2 + k 2 ) = ^-G(±d(a)pl + a 3 V{<j))) (17) 

and from the equations of motion generated by the constraint the effective Raychaudhuri equation 

d 8ttG ( , . 9 / -i dlog(a 3 d(a))\ T „. A , 10N 
a = [ a d(a)P * V - da ) ~ VW ) ' (l8) 

In the latter, the modification leads to an entirely new term. 



3 Phenomenology 

These equations, all resulting from a single modification in the matter Hamiltonian implied by effects of 
the basic quantum representation, are the starting point for phenomenology of loop cosmology based on 
effective equations. Additional corrections, related to quantum fluctuations and additional quantum degrees 
of freedom, have not yet been studied systematically but are expected to be important only on very small 
scales. Modifications in the matter Hamiltonian, on the other hand, are non-perturbative and can be shifted 
into regimes of larger scales just by choosing a large value for j. This allows one to study its implications 
in isolation from other correction terms, even though one does not expect j to be very large. Additional 
corrections |48 47| and similar parameters |5l | also arise for the gravitational part of the constraint and 
become relevant when the matter density is large (of Planck size). 

The main modification can then easily be interpreted intuitively by viewing the Friedmann equation as 
the energy equation of a classical mechanics system with a potential determined by the matter Hamiltonian. 
Classically, the matter Hamiltonian is usually decreasing as a function of a as a consequence of the kinetic 
term. This is a consequence of the fact that classical gravity is always attractive. When a -3 is replaced by 
d(a) as in Fig.^ however, the slope of the potential is flipped on small scales and it becomes increasing. 
Thus, the direction of the force changes and (quantum) gravity on small scales receives a repulsive contri- 
bution. With this picture, one can easily imagine that characteristic effects can arise such as the prevention 
of collapse into a singularity by a bounce or the acceleration of expanding evolution to an inflationary era. 



3.1 Bounces 

For a bounce to be realized, we need to find a time where a = and a > 0. The first condition can 
be checked with the effective Friedmann equation (Il7> where, for d = to be possible, the positive 
kinetic term must be compensated by either a positive curvature term k = l I53I I54I or a negative scalar 
potential V(<p) < |55|. The second condition at solutions for a — is then controlled by the effective 
Raychaudhuri equation ( Il8> which for classical solutions gives only negative d, i.e. recollapse points. 

When a turning point falls in the modified regime, on the other hand, the additional term in the effective 
Raychaudhuri equation as well as modified matter behavior through the effective Klein-Gordon equation 
d!6i can change the picture and imply bounces. This is not realized in all cases, but happens generically 
and without the need for special potentials in contrast to the classical situation. Numerical solutions for 
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two examples are given in Fig. [5] These effective bounces can be seen as a consequence of a repulsive 
gravitational force, or equivalently of negative pressure 



P = - 



dH 
W 



1 dH 

3a 2 do 



< 



(19) 



which changes sign on scales where the matter Hamiltonian is increasing. 
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Fig. 3 Numerical solutions (solid lines) for bounces in models with positive spatial curvature (left) and a negative 
potential (right), respectively, compared to the classically singular solutions (dashed). 



The modification in the matter Hamiltonian alone does not give rise to bounces in all cases, such as a 
flat model with positive potential. A universe then still collapses to small sizes closer and closer to the deep 
quantum regime. Eventually, its dynamics has to be described by the basic difference equation which is 
non-singular, but at small sizes also additional corrections become important in effective equations. Those 
corrections are indeed directly related to the underlying discreteness or quantum fluctuations of geometry, 
and appear on very small scales. This, then, provides bounces more generically 1 50 56, 49ll21l . 



3.2 Inflation 

Negative pressure is also required for inflation where a > can now follow without special conditions for 
the potential or the initial values of an inflaton field. In fact, inflation in the general sense of accelerated 
expansion now happens generically and even without any potential at all [52 1 as illustrated in Fig.|4] It 
is enough that the kinetic term in the effective Friedmann equation becomes increasing as a function of a 
which is always realized on small scales. 

The precise manner depends on the ambiguity parameter I for which we obtain d(a) ~ a?^- l ~ 1 ' on small 
scales which is always increasing due to < / < 1. This form determines, for a vanishing potential, the 
type of inflation which is super-inflationary (with equation of state parameter w = — 1/(1 — /) < —1). For 
more realistic models with a potential, however, the behavior is driven very close to exponential inflation 
thanks to the a-dependence of the potential term 1571 . From the behavior of matter in such an inflating 
background one can then derive potentially observable effects. 

We first look at the kinetic term driven inflation which happens for any matter field and thus can elimi- 
nate the need for an inflaton. During accelerated expansion in this regime, structure can be generated with 
less fine tuning than in inflaton models. As preliminary calculations indicate [58 1, the resulting spectrum is 
nearly scale invariant 1 57 1 , and quantum effects can also be used to give arguments, based on 1 59 60 1 , for a 
small amplitude in agreement with observations. Details of the spectrum such as the running of the spectral 
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Fig. 4 Numerical solution with inflationary behavior (right) in the regime of modified effective densities 
(left). Dashed lines indicate the range of sizable modifications. 

index do depend on ambiguity parameters such that they may be restricted by forthcoming detailed obser- 
vations. In contrast to single field infiaton models, the spectral index would be a little larger than one as a 
consequence of super-inflation, i.e. the spectrum is slightly blue, which can be a characteristic signature. 

If this phase alone would have to be responsible for the generation of all structure, an extremely large 
value for the parameter j would be required. Moreover, the spectrum would then be blue on all scales 
which is ruled out by observations. As we will see now, however, the first inflationary phase will always 
be followed by slow-roll phases because matter (or infiaton) fields are driven away from their potential 
minima while the modified Klein-Gordon equation is active. Subsequent phases can thus serve to make 
the universe large enough and provide a sufficient amount of e-foldings, while visible structure can have 
been generated in the first, quantum phase. Attractive features of this class of scenarios are that potentials 
for matter fields driving late-time slow-roll phases do not need to be special and that visible structure would 
have been generated by the quantum phase, thus giving potentially observable quantum signatures. 

3.3 Scalar dynamics 

In the effective Klein-Gordon equation dl6l the classical friction term of an expanding universe, used for 
slow -roll behavior, changes sign and turns into an antifriction term on small scales 1521 . Matter fields in 
this regime are then excited and can move up their potential walls even if they start close to minima at 
small momenta. After antifriction subsides, they will continue to roll up the walls but be slowed down by 
the then active friction. Eventually, they turn around and roll down the potential slowly. In this manner, 
additional phases of inflation are generated. The whole history is shown in Fig.[5]with a rapid initial phase 
containing the push of tj> up its potential and the loop inflationary phase together with the later slow-roll 
stage. A solution for the scale factor showing both inflationary phases in the same plot is given in Fig. [6] 

If <f> is an infiaton such as that of chaotic inflation with a quadratic potential, initial conditions are 
provided for a long phase of slow-roll inflation. But the slow-roll conditions are not satisfied in all stages, 
in particular not around the turning point of the infiaton. Here, <f> vanishes or is very small such that the 
slow -roll condition cf) <C H<j> cannot be satisfied. If the second phase is responsible for structure formation, 
structure on large scales, generated in early stages of the slow-roll phase, will differ from usual scenarios. 
This can in particular contribute to a suppression of power on large scales |61 1. If the last inflationary 
phase takes too long these effects are not visible, but one can estimate the duration from the infiaton initial 
values one typically gets through the antifriction mechanism. As it turns out, one often obtains observable 
effects, i.e. there is a sufficient amount of inflation but not too much for washing away quantum gravity 
effects |62|. This aspect is even enhanced if one looks at the evolution at the level of difference equations 
(63). 
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Fig. 6 Numerical solution showing the first inflationary phase of loop cosmology and the beginning of the 
ensuing slow-roll phase after the scalar is pushed up its potential. 



3.4 Combinations of different phases 

So far we have considered individual bounces or inflationary phases. Depending on details of the matter 
system such as potentials, many combinations are possible and often of interest for model building. If 
one combines several bounces in models with a classical recollapse, oscillatory models result [64]. The 
behavior can then gradually change from cycle to cycle and eventually, after small changes have added up, 
result in a qualitative change in the universe behavior. 
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This is realized, for instance, in a new version of the emergent universe which was originally devised 
as a non-singular inflationary model with positive spatial curvature |65|. The singularity is avoided by 
starting the model close to a static Einstein space in the infinite past which, with a suitable potential, 
develops into an inflationary phase. The initial state is, however, very special because the static Einstein 
space is unstable. This changes if effective equations of loop cosmology are used: there are new static 
solutions on small scales which, in contrast to the classical ones, are stable [66|. Starting close to those 
solutions will lead to a series of cycles of a small universe which can gradually change due to the motion 
of matter fields in their potential. 

In this manner, one obtains a whole history for a universe evolving through different phases of con- 
traction and expansion in a cyclic but not necessarily periodic manner. However, the cycles are usually 
short and do not automatically give rise to a large universe. On the other hand, matter fields move in their 
potential during the evolution, and close to bounce points antifriction effects can help to bring matter fields 
far up their potentials or over potential barriers. If a suitable part of the potential is encountered during 
this process, the conditions can be right for the start of a slow-roll regime of sufficient duration for a large 
universe to emerge |67|. In fact, this can be obtained in the emergent universe scenario from a simple 
initial state evolving, after many cycles, to an inflationary phase for the structure formation we need for 
our universe |66 68 1. A numerical solution for this scenario is illustrated in Fig.0 The characteristic 
feature is that it is based in an essential manner on closed spatial slices of positive curvature which, if the 
inflationary phase is not too long, may be detectable in the near future |66|. 
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Fig. 7 Stroboscopic density plot of solutions in the classical phase space (arbitrary units) with initial 
cyclic behavior and an eventual inflationary phase. To illustrate the time behavior, Gaussians have been 
added peaked at discrete coordinate time intervals. 

When different matter sources are present, such as different fields or a fluid in addition to a scalar field, 
other possibilities arise. There can still be fixed points which allow cyclic behavior around them, but the 
position can now depend on the field values. In particular, cyclic motion of the fixed points themselves is 
possible which implies that the evolution is double-cyclic with small cycles around a fixed point superposed 
to the cyclic motion of the fixed point 1 69 1 . There are several new possibilities which are being investigated 
in dynamical system approaches. In this context it is in particular of interest when a graceful entrance 
into an inflationary phase can arise in the presence of different matter sources, i.e. how generically initial 
conditions become right to start a sufficiently long slow-roll phase. 
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Cyclic behavior is also studied often in the context of brane collisions. While loop quantum gravity 
and cosmology are difficult to formulate in this higher-dimensional setting, 1 one can simply take potentials 
motivated from brane scenarios 11711 and study implications with loop modifications. The interpretation of 
the scalar is then as a radion field, i.e. the distance between two branes. A common characteristic property 
is the possibility of negative potentials which, as we have seen earlier, allow bounces in loop cosmology. 

Loop cosmology then provides a non-singular bounce for such a brane model, and one can check how 
easy it is to realize properties assumed in other models where a mechanism for singularity removal was 
not known. In some of those scenarios, for instance, it was assumed that the scale factor bounces and 
simultaneously the scalar field turns around 11721 . Only then does one really have a model where the branes 
first approach each other and then bounce off. This turns out to be impossible in loop cosmology based 
on dl7> where bounces are realized with a negative potential but, as illustrated in Fig. [8] the scalar cannot 
change direction |55|. The only possibility for <fi to turn around is if it encounters a region of positive 
potential, but then one automatically obtains an inflationary phase and the scenario is not different from 
those described before. Similar results can be obtained when bouncing solutions realized with higher 
curvature corrections are used 1731 . 




Fig. 8 Bouncing numerical solution for the scale factor (left) in a brane potential with a monotonic decrease of 
the scalar (right). During bounces, the scalar is only slowed down due to antifriction in the effective Klein-Gordon 
equation, but it cannot turn around. A turnaround can be realized only if the scalar encounters positive values of its 
potential, which generically leads to inflation (dashed). 



4 Conclusions 

Through its background independent and non-perturbative quantization scheme, loop quantum gravity 
and cosmology are well-equipped for extreme physical situations as they are realized close to classical 
singularities. With symmetric models, those situations can now be dealt with in many cases including 
cosmology and black holes. Many examples exist for non-singular models by the same general mechanism, 
and also conceptual problems can be solved. In addition, many phenomenological applications arise from 
a few basic effects. 

These effects have, in all cases, been known first from mathematical considerations and then transferred 
to and evaluated explicitly in models. The discussed effects then resulted automatically, rather than being 
looked for with particular applications in mind. Moreover, a few effects going back to properties of the ba- 
sic representation dictated by background independence suffice to cover a plethora of physical applications 
in different areas. 



There is, however, an interesting duality to brane-world models if conditions for the shape of d(a) are relaxed 70]. 
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This fact is encouraging for the viability of the whole framework, but it is still important to understand 
the relation between models and the full theory as completely as possible. Some properties can directly be 
related, and others have to be derived after analogous constructions. There are no known contradictions, 
but at the dynamical level also no proof of, e.g., singularity removal without assuming symmetries. It is, 
however, clear already that the presence of local physical degrees of freedom is by itself no obstacle to 
non-singular quantum evolution. 

Also models, in particular inhomogeneous ones, have to be developed in more detail. At the fun- 
damental level one then has to deal with many coupled partial difference equations, for which even the 
formulation of a well-posed initial or boundary value problem can be difficult. For solutions one will have 
to refer to new numerical techniques which are being developed in investigations of numerical quantum 
gravity. Effective equations would also help considerably in understanding inhomogeneous situations, but 
their derivation is much more complicated than in homogeneous models and so far not completed. 

At the same time, properties of the full theory are being understood better. There is thus an approach 
from two sides, by weakening symmetries in models to get closer to the full setting and by starting to 
understand full configurations which can be argued to be close to states considered in a symmetric model. 
A third direction will, in the future, be provided by observations and their relation to phenomenological 
results. In this way, confidence in physical effects derived with loop methods will be strengthened and can 
eventually be compared with observations. 
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